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ABSTRACT 

^ We present a detailed study of the emission from a z — 2.35 galaxy that causes damped 

^ Lyman-a absorption in the spectrum of the background QSO, SDSS J 2222—0946. We present 

the results of extensive analyses of the stellar continuum covering the rest frame optical- 
1-1 UV regime based on broad-band HST imaging, and of spectroscopy of the strong emission 

lines: Lya, [Oil], [Olll], [Nil], and the Ha and H/3 lines from VLT/X- Shooter. We compare 

the metallicity from the absorption lines in the QSO spectrum with the oxygen abundance 

I inferred from the strong-line methods (R23 and N2). The two emission-line methods yield 

£T} consistent results: [O/H] = —0.30 ± 0.13. Based on the absorption lines in the QSO spectrum 

{vq a metallicity of— 0.52±0.05is inferred at an impact parameter of 6.3 kpc from the centre 

t^J- of the galaxy. The star formation rates of the galaxy from the UV continuum and Ha line 

can be reconciled assuming an amount of reddening of E(B — V) = 0.06 ± 0.01, giving 

an inferred SFR of 13 ± 1 M yr _1 (Chabrier IMF). From the HST imaging, the galaxy 

^£) associated with the absorption is found to be a compact (r e =1.12 kpc) object with a disc- 

tH like, elongated (axis ratio 0.17) structure indicating that the galaxy is seen close to edge on. 

r* Moreover, the absorbing gas is located almost perpendicularly above the disc of the galaxy 

. J^J suggesting that the gas causing the absorption is not co-rotating with the disc. We investigate 

K> the stellar and dynamical masses from SED-fitting and emission-line widths, respectively, and 

^ find consistent results of 2 x 10 9 M . We suggest that the galaxy is a young proto-disc with 

evidence for a galactic outflow of enriched gas which cools down at a large distance. This 
galaxy hints at how star-forming galaxies may be linked to the elusive population of damped 
Lyman-a absorbers. 

Key words: galaxies: formation - galaxies: high-redshift - galaxies: ISM - quasars: absorp- 
tion lines - quasars: individual: SDSS J 22:22:56.1—09:46:36.2 - cosmology: observations 



1 INTRODUCTION 

Mapping the structure, properties and chemical enrichment of 
* Based on observations collected at the European Organisation for As- galaxies over cosmic history is a major goal of contemporary as- 
tronomical Research in the Southern Hemisphere, Chile, under program trophysics . T n me local Universe important advances can be made 
087.A-0085(A). Based on observations made with the NASA/ESA Hubble , /, . ., ^ 1V . , J c + • u. o 1 • u 

rT , 11 . 11f , m t „ . T . , . , by studying the age-metallicity relation 01 stars in the Solar neish- 
Space Telescope, obtained at the Space Telescope Science Institute, which ~ Trr~\ — r ■ ■■ ■ ■ 

is operated by the Association of Universities for Research in Astronomy, ^ '°''l — — ° ' ' - — J — ; — 1 



Inc., under NASA contract NAS 5-26555. These observations are associ- g rou P dwarf galaxies JFrebeletal. , 2007. 2010). At higher redshifts 



ated with program 12553. than z ~ 2, a powerful method for studying chemical evolution 

f E-mail: krogager@dark-cosmology.dk is spectroscopy of Damped Lya Absorbers (DLAs) detected either 



2 J.-K. Krogager et al. 



towards background QSOs ( [Pettini et al.|1990[|Wo?fe et al.|2 005, 
M0ll er et al.|2013|) or against Gamma-ray Burst (GRB) afterglow 
light jFynbo et al.[[2006l |Savaglio|[2006l |Prochaska efaTpOOT] ). 
In those studies, abundances are determined from the gas phase as 
probed by the H I and metal absorption lines detected against the 
light of the background QSOs and GRB afterglows, respectively. 

Mainly at lower redshifts, but recently also at z > 2, H II- 
region abundances are determined using the relative strengths of 
strong emission lines (e.g., Kewley & Ellison 2008, Shapley 2011 , 
and references therein). Currently, there have only been very few 
cases where both methods have been applied to the same object. 
The first example is the case of the DLA towards SBS 1543+593 
(Bow en et al.|2 005) where the two methods yielded consistent re- 
sults. Since then a handful of other cases have been studied (see the 
compilation in Peroux et al. 20 j^|. As seen in Fig. 8 of Peroux etal.| 
( |2012| ), the absorption line measurements generally probe regions 
at larger galactocentric distances than the emission-line based mea- 
surements. The absorption-line measurements on average indicate 
lower metallicities than the emission-line based measurements pos- 
sibly reflecting the early setup of metallicity gradients ( O'Rourke 
|et al.|20lT) . However, it is important to stress that different recipes 
to derive oxygen abundances using strong emission-line fluxes re- 
veal very inconsistent results (Kewley & Ellison 2008 ). It is there- 
fore of interest to expand the sample of sources where strong-line 
based abundances can be independently tested using other methods 
(see also |Pettini|2006l|Kudritzki et al.|20T2| ). Also the use of dif- 
ferent elements, e.g., Fe, Zn or Si, to infer absorption metallicities 
may introduce systematic offsets when comparing galaxies within 
a heterogeneous sample. 

Combining the two complementary methods of studying the 
metal enrichment in galaxies provides important hints to under- 
standing how galaxies turn their gas into stars, as the absorption 
lines directly probe the cold gas, and the properties of the star 
forming region can be probed directly from the emission lines, and 
sometimes also from the continuum. However, linking the absorp- 
tion characteristics of DLAs to the emission characteristics of the 
galaxies causing the absorption has been a great challenge at high 
redshift, due to the faint nature of DLA galaxies and due to their 
proximity to a very bright QSO. So far, only a few DLA systems 
with confirmed emission counterparts have been established at red- 
shifts around and higher than two (Krogager e t al.|2012| . 

Here we present an analysis of the DLA at z = 2.35 seen in 
the spectrum of z = 2.93 QSO 2222—0946. The galaxy counter- 
part of the absorber was detected by Fynbo et al. (2010 ) and here we 
present new, deeper spectroscopic data allowing us to measure the 
metallicity directly from nebular [Oil], [Olll], and B aimer emis- 
sion lines. Furthermore, we combine the spectroscopic data with 
imaging from Wide Field Camera 3 (WFC3) onboard the Hub- 
ble Space Telescope (HST). The high resolution images from HST 
makes it possible to detect the continuum emission from the galaxy 
directly, allowing us to characterise the properties and structure of 
the absorbing galaxy. This enables us to start bridging the gap be- 
tween the population of DLA galaxies and star forming galaxies. 

The paper is organised as follows. In Sect. 2, we present 
our spectroscopic data from X-Shooter and our imaging data from 
HST/WFC3. In Sect. 3, we derive the oxygen abundance from neb- 
ular emission lines, derive the gas phase metal abundances, char- 
acterise the morphology from the imaging data, and finally we 
obtain fluxes of the galaxy used for spectral energy distribution 
(SED) fitting. In Sect. 4, we compare the absorption and emission 
properties, and discuss and outline the implications of our work. 



Throughout this paper, we assume a standard ACDM cosmology 
with Ho = 71 kms^Mpc" 1 , Ov = 0.7 and Q M = 0.3. 



2 OBSERVATIONS AND DATA REDUCTION 

2.1 X-Shooter spectroscopy 

QSO 2222—0946 was observed with the X-Shooter spectrograph 
( Vernet et al. 201 1 ) mounted on ESO's Very Large Telescope, Unit 
Telescope 2, for a total of 10 hours and 47 minutes during the 
nights of 201 1 May 14, 201 1 August 4, 29, and 30. The integrations 
were split into ten different exposures each made up of four sub- 
exposures dithered along the slit. For all observations slit- widths of 
0.8, 0.7 and 0.9 arcsec for were used for UVB, VIS and NIR, re- 
spectively, and we used the slow readout with 1x2 binning, i.e., 2 
times binning in the spatial direction. The individual spectra have 
been processed using the official ESO pipeline (Modigliani et al.| 
|2010| ) for nodded exposures version (2.0.0). The spectra were in- 
dividually flux-calibrated using the standard star observed at the 
beginning of the night for each integration. All spectra were taken 
during very good conditions with clear sky and good seeing. The 
error on the flux calibration is on the order of 5%, as determined 
from the robustness of the calibration when using different stan- 
dard stars. We note that Fyn bo et al.| ( |2010| ) scale their spectrum 
of this target to the one available from the Sloan Digital Sky Sur- 
vey. Given our very robust flux calibration we have chosen to trust 
the pipeline product. Also, comparison of spectra observed in very 
different epochs, as is the case here, is not straightforward given 
the random variability of QSOs. We encountered issues with a too 
low flux-level of the VIS spectra for unknown reasons. For this 
reason we scaled the VIS spectra to the overlapping regions of the 
well-determined UVB and NIR spectra. With a simple multiplica- 
tive factor we were able to match both UVB and NIR at the same 
time. The scaling does not affect any of the extracted emission lines 
as they are all in the NIR spectra, except for Lya in the UVB. 

The final 2D- spectra were then spatially aligned and added 
by error- weighting, and we correct for slit-loss by calculating how 
much light gets dispersed outside the slit at the given seeing in each 
exposure before co-adding all the observations. We have not cor- 
rected for telluric absorption, since those regions are not crucial 
for our analysis. The effective seeing in the final flux and wave- 
length calibrated 2D- spectrum is 0'/7 (inferred from the width of 
the trace at 715 nm) and we measure the resolving power in the 
spectrum from telluric lines to be 1 1,000 and 7,000 for the VIS and 
NIR arms, respectively. As there are no direct features in the UVB 
arm to indicate the resolving power we assume the resolving power 
in the UVB arm to be the expected 6,200 from the specifications 
of the X-Shooter manual given a slit width of 0V8, as the average 
seeing in UVB was / /77 (measured at 480 nm). 

In Figure[T]the spectrum blue- ward of the QSO hya emission 
line is shown, demonstrating the quality of the combined spectrum; 
note the significant flux below the Lyman limit at z = 2.87. 



2.2 HST I WFC3 imaging 

The field was observed with the Wide Field Camera 3 on 201 1 Nov 
10 (with the UVIS detector in the F606W filter) and on 2012 Sep 14 
and 15 (with the IR detector in the F105W and F160W filters). The 
roll-angle of the telescope was set such that the galaxy counterpart 
of the DLA fell between the diffraction spikes of the Point Spread 
Function (PSF). The two observations with IR detector were taken 
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using the WFC3-IR-DITHER-B0X-MIN pattern providing an op- 
timal 4-point sampling of the PSF. The UVIS observation was taken 
using the WFC3 -UVIS -DITHER-BOX pattern. 

The field was observed using a 4-point sub-pixel dither pat- 
tern, which allowed us to regain more spatial information. We have 
used the software package multidrizzle[jto align and com- 
bine the images. By shifting and combining the images taken with 
sub-pixel offsets one achieves a better sampling of the PSF, which 
in the case of the NIR observations is quite crucial as the PSF is 
poorly sampled in the native 0V13 px _1 images. Furthermore, the 
drizzle- algorithm allows us to reduce the pixel size when combin- 
ing the images providing sharper and better sampled images. For 
the combination in this work we have set the parameter pixf rac 
to 0.7 and used a final pixel scale of 0V06 px _1 for NIR and 
/ /024 px _1 for UVIS. For a detailed description of the parame- 
ters in the software we refer the user to the mult idr i z z le 1 user 
manual. 



fit parameters. The Lya line is shown in Fig. [3] Due to the asym- 
metric shape and the complex nature of the resonant line, we deter- 
mine the line flux simply by integrating the observed line profile as 
opposed to fitting the line. 

3.3 Balmer Decrement 

The ratio of line the Ha and H/3 line fluxes provides us with in- 
formation about the dust extinction in the system as we know what 
this ratio should be intrinsically, given the physical conditions of 
the emitting region. Assuming case B recombination with an elec- 
tron temperature of T e — 10 4 K and density of n e — 10 2 cm -3 , 
which are standard assumptions in the literature for star forming 
regions, the line ratio has a value of 2.86. Requiring that our mea- 
sured line ratio, after reddening, and the intrinsic ratio be the same, 
we estimate the reddening: 



E(B - V) 



2.5 



fc(H/3) - k(Ra 



log 



(Ha/H/3) obs 
(Ha/H/3) 



3 RESULTS 

3.1 Spectral PSF subtraction 

In order to detect the faint emission lines from the foreground 
Damped Lyman-a Absorbing galaxy we needed to subtract the 
QSO continuum. We did this by modelling the spectral trace and 
the spectral point spread function (SPSF) as a function of wave- 
length. The trace position and amplitude were fitted as a function of 
wavelength in a small (84 A) region around each emission line from 
the galaxy, e.g., [O II], [O III], and [N n]. In the same region around 
the line we estimated the SPSF by averaging the observed spatial 
profile along the spectral direction. The modelled 2D- spectrum of 
the QSO was then subtracted from the observed spectrum, and a 
ID-spectrum for each line was extracted with the optimal extrac- 
tion algorithm of Home (1986). All extracted emission lines are 
shown in Figure [2] 



3.2 Extracting the Emission Lines 

We were able to detect emission from Lya, Ha, H/3, the 
[O II] AA 3727, 3729 doublet, the two [O III] lines at 4959 A and 
5007 A, and [Nil] A 6583. The line profiles for all the extracted 
lines were fitted with Gaussian profiles to measure the flux in the 
line. The parameters of the Gaussian were allowed to vary in all 
fits except for [O n] and [N n] . The line- width and redshift of the 
two [O n] components were tied in the fit, because there is a sky 
line that falls right between the two components. For the very faint 
[N n] -line, we re-binned the spectrum by a factor of two and fixed 
the redshift and line-width to the values of the nearby Ha line. 
The errors on the fluxes were estimated by varying the line profile 
within the error of each fit parameter 1,000 times. The uncertainty 
was then determined from the 1 a width of the resulting distribution 
of fluxes. The fluxes and errors are listed in Table[T] 

All the extracted lines with Gaussian line fits are shown in 
Fig.[2]along with the best fitting gaussian profile. The grey shaded 
areas indicate sky lines or telluric features. These regions were in- 
cluded in the fits to give the optimal estimate of uncertainties on the 



1 multidrizzle is a product of the Space Telescope Science Institute, 
which is operated by AURA for NASA. 



where k denotes the reddening law evaluated at the given wave- 
length and (Ha/H/5)o indicates the intrinsic line ratio. We use the 
extinction law from Calzetti et al. ('2000) to quantify the extinction 
while adopting a Rv value of 4.05. From the measured ratio in our 
spectrum of (Ha/H/3) bs = 3.03 ± 0.34 we obtain the corre- 
sponding extinction of E(B — V) = 0.05±0.10. In this calculation 
we have not taken into account effects from differential slit-loss, 
which arise from the wavelength dependence of the seeing and at- 
mospheric dispersion. The slit loss is greater at lower wavelengths, 
i.e., for the H/3 line in the case of the Balmer line ratio. However, 
the effect in the K- and if -band, where we extract the Balmer lines, 
is minor (~ 2%) compared to the uncertainty introduced by the sky 
lines, especially at the position of the H/3-line. 



3.4 Metallicity 

We have inferred the metallicity of the system using three indepen- 
dent methods; two measures of the emission-line metallicity using 
the strong line ratios R23 and N2, respectively, and one measure of 
the absorption-line metallicity from Voigt-profile fitting. In the next 
sections we present each determination in detail. 



3.4.1 R23 calibration 

We determined the metallicity 12 ± log(0/H) by use of the 
strong line diagnostic, R23 ( |Pagel et aL]|1979| >, defined as the ra- 
tio ([On] AA 3727, 3729 + [Om] A 4959 + [Om] A 5007) / H/3. 
One complication of using this diagnostic is that the metallic- 
ity is double-valued for a given value of R23 (the so-called up- 
per and lower branches). However, due to our relatively high 
value of R23 we are in the region of the diagram where the two 



branches are close to one another. The calibration by Kobulnicky 
& Kewley (2004) that we use depends on the ionisation param- 
eter, q, which can be determined by using the line ratio O32 = 
[O III] A5007 / [0 11] A3727, but requires prior knowledge of the 
metallicity. Given the two measured quantities, R23 and O32, we 
can solve for both metallicity and q at the same time by iterating 
the two expressions given an initial guess of metallicity. 

From our spectrum we measure log(R23) = 1.00 ± 0.05 and 
log(032) = 0.55 ±0.05. These line ratios are corrected for redden- 
ing as indicated by the Balmer decrement. However, in Sect. |3.8| we 
infer E{B — V) — 0.06 ±0.01, which is the value that we have 
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Table 1. Measured emission line fluxes 



Transition 


Wavelength^ 1 ) 


fW 2 ) 


fwhm( 3 ) 


z 


Lya 


1215.67 


14.3±0.3 


- 


- 


[On] 


3726.03,3728.82 


2.9±0.3 


121 ±8 


2.3536 


H/3 


4861.33 


1.9±0.2 


151 ± 15 


2.3537 


[Om] 


4958.92 


4.1±0.3 


129 ±5 


2.3537 


[Om] 


5006.84 


11.6±0.6 


110 ±2 


2.3537 


Ua 


6562.80 


5.7±0.3 


124 ±4 


2.3537 


[Nil] 


6583.41 


0.6±0.2 


- 


- 



t 1 ) Transition rest frame wavelength in A. 

( 2 ) Flux in units of 10 — 17 erg s — 1 cm -2 , before reddening correction. 

( 3 ) Line width at FWHM in units of km s _1 corrected for the instrumental 
resolution of 45 km s — 1 . 

adopted for the correction here. We perform the iterative calcula- 
tion of metallicity 500 times while varying the two line-ratios (R23 
and O32) within their respective errors. The metallicity and uncer- 
tainty is then determined as the median and standard deviation of 
the resulting distribution This in turn yields an inferred metallic- 
ity of 12 + log(0/H) = 8.41 ± 0.19, see Fig.g We furthermore 
note that the result does not depend on the initial guess, and that 
the model uncertainty introduced by the error on the ionisation pa- 
rameter, q, is negligible. The uncertainty of 0.19 dex includes the 
contribution from the scatter in the R23 calibration which we con- 
servatively set to 0.15 dex. 



3.4.2 N2 calibration 

Our detection of [N 11] enables us to infer the metallicity indepen- 
dently from the N2 index ( [N 11] / Ha ) as this line ratio has been 
found to correlate with metallicity (Pettini & Pagel 2004 ). The cal- 
ibration from these authors give 12 + log(0/H) = 8.90 + 0.57 • 
log(N2) with a scatter of ±0.18 dex. The small wavelength differ- 
ence of Ha and [N 11] means that we do not need to worry about 
reddening effects in the N2 ratio. The measured line fluxes yield a 
line ratio of N2 = 0.11 db 0.03, giving an inferred metallicity of 
12 + log(Q/H) = 8.36 ±0.19. 

Peroux et al.|(2012) infer an upper limit of 12 + log(0/H) < 8.2 
from their non-detection of [Nil]. Our measurement is slightly 
higher, but still consistent with theirs within the uncertainties. 

Since the two measures of the emission metallicity, R23 and 
N2, are independent we can combine both measures to obtain a 
more precise estimate of 12 + log(0/H) = 8.39 ± 0.13. Adopting 
the Solar abundance presented by Asplu nd et al. 1(2009) this gives 
us [0/H] = -0.30 ±0.13. 



3.5 Absorption-line Abundances 

In the following section we present detailed measurements of the 
absorbing abundances detected in the DLA system. |Fynbo et al.| 
( |2010| ) present an analysis of the metal absorption-line system of 
this DLA using a three-component Voigt profile fit. They derive 
a fairly high metallicity from [Zn/H] = -0.46 ± 0.07, [Si/H] = 
-0.51 ± 0.06, [S/H] = -0.6 ± 0.1, [Mn/H] = -1.23 ± 0.06, 
[Ni/H] = -0.85 ± 0.06, and [Fe/H] = -0.99 ± 0.06. Since, 
in this work, our data have much better resolution (recall values 
in Section [TT] compared to R = 4700, 6700, 4400 in the Fynbo 
et al. study for the UVB,VIS and NIR arms, respectively) and 
much higher Signal-to-Noise Ratio (SNR^lOO, compared to their 
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Figure 2. Emission lines extracted after the subtraction of the QSO contin- 
uum in the 2D spectrum. Each panel shows the 2D spectrum where the hor- 
izontal residuals from the QSO subtraction is seen (top) and the extracted 
ID spectrum of the line (bottom). The wavelengths are given in /im and all 
flux-density units are 10 — 17 erg s — 1 cm -2 A . The blue line in each 
panel indicates the best fitting Gaussian to the line profile. The grey filled 
areas indicate skylines or telluric absorption features. 



SNR^50), we can alleviate the effect of hidden saturation, thereby 
allowing us to validate the previous measurements. 

We have focused on Zn II, Si II, and Fe II for the purpose of 
determining the metallicity. In order to properly decompose the ab- 
sorption profiles in different velocity components we have selected 
weak, but well-defined lines. We found that a five-component fit 
provides a satisfactory description of the absorption lines; how- 
ever, we had to include the Mg I AA 2026,2852 lines in order to fit 
the Zn II A 2026 line. All the lines were fitted simultaneously and 
the results are summarised in Table |2] We derive the total column 
densities from the fit by adding the column densities for all com- 
ponents of each individual line. As a result we obtain the following 
metallicities: [Zn/H] = -0.51 ± 0.06, [Si/H] = -0.53 ± 0.05, 
and [Fe/H] = — 1.02 ± 0.05, assuming Solar abundances from |As-| 
plund et al. ('2009). We do not consider the abundance of Mg as 
most of this element will be in form of Mg II, and the fit to the 
Mg I A 2852 line is influenced mildly by telluric absorption. The 
inferred metallicities are in excellent agreement with those deter- 
mined by Fynbo et al. (2010). We have used the column density of 
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Figure 1. The final ID spectrum from the UVB arm of the z = 2.9 Q 2222—0946. In the left part of the spectrum, the Lyman limit from an absorber at redshift 
z = 2.87 is clearly seen and the damped Lyman-alpha absorption line is identified in the middle of the figure by the vertical dashed line. The Lyman-alpha 
emission line at the QSO redshift is seen in the right part of the figure. 




Table 2. Ionic column densities for individual absorption components of 
the z = 2.3533 DLA. The velocity indicated at each component shows the 
relative velocity offset with respect to the adopted systemic redshift. 



4070 4080 

Wavelength [A] 



Figure 3. The Lycn emission line extracted from the DLA trough in the X- 
shooter spectrum. The top panel shows the observed 2D spectrum and the 
panel below shows the extracted ID spectrum, and la uncertainty (grey 
shaded area). The dotted vertical line shows the systemic line centre. The 
common asymmetric line shape of Lycn is seen in this case, where the blue- 
shifted part of the line is heavily absorbed. On top of the observed spectrum 
the best fitting model of the emission profile is shown (red line) along with 
the 1 a confidence interval (red shaded area), see Sect.|3.7|for details. 



log(Nm) = 20.65 ± 0.05 from the Lya absorption line of. The 
uncertainty of our metal column densities, which we estimate to 
be 0.05 dex, is dominated by the uncertainty on the normalisation 
of the QSO spectrum. The inferred metallicity places this DLAs 
among the < 5% most metal rich DLAs known today (e.g., |Noter^| 
|daeme et al.p008] ). 
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Figure 4. Plot of the R23 strong line ratio using the calibration by |Kobul^] 
nicky & Kewley ( 2004 ) (two solid lines) with model uncertainty of 0.15 dex 
indicated by the light grey shaded area around the curves. The orange 
filled area between the two solid lines indicates the constraint from our 
measurement of R23 . The dashed line indicates the inferred metallicity of 
12 + log(0/H) = 8.41 and the red shaded area around it indicates the 
lcr uncertainty of 0.19 dex. The dot-dashed line shows the Solar oxygen 
abundance of 12 + log(O/H) = 8.69 jAsplund et al.|2009) . 
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Figure 5. Results of Voigt-profile fitting to the low-ionization metal ab- 
sorption lines from the z = 2.353 DLA towards Q 2222-0946. The 
Mg I A 2852 line is influenced mildly by telluric absorption lines. 



3.6 Structural Fitting of the HST images 

We started out by modelling the PSF in each image using the 
software TinyTim to create sub-sampled simulated Point Spread 
Functions for each of the individual frames each of the positions 
in the dither pattern. We preferred this approach over using stellar 
sources as the TinyTim PSFs are more sensitive to the outer re- 



gions of the PSF, where the stellar sources are dominated by noise 
in the sky background. 

After resampling the modelled PSFs to the native pixel size and 
convolving with the filter- specific Charge Diffusion Kernel, these 
"raw" PSF images were drizzled together in the same way as the 
actual data to replicate the effect of multidrizzle on the PSF 
shape. We were not able to simulate the PSFs for the infrared im- 
ages properly because the images were mildly saturated. We there- 
fore chose to use stellar PSFs generated by median combination of 
stars in the field. Since there are only a handful isolated stars in the 
field of view with signal-to-noise ratios similar to that of the QSO, 
we used a few stars with slightly lower SNR in order to reduce the 
background noise in the PSF. 

We used the software Galfit ( |Peng et al.|2002| ) to subtract 
the QSO and to characterise the absorbing galaxy. Using the mod- 
elled PSF for the F606W image we first subtracted the QSO by 
simply modelling it as a point source while modelling a constant 
sky background. We then located the nearby galaxy in the residu- 
als. Hereafter, we re-did the fit, this time simultaneously fitting the 
QSO, the background, and the galaxy using a Sersic surface bright- 
ness profile: 



E(r) = E e exp 



l/n 



1 



where H e is the surface brightness at the effective radius, r e , de- 
fined as the radius enclosing half the flux. The Sersic index, n, 
determines the concentration of the profile, with high n profiles 
having steeper inner slopes and larger extended wings. The oppo- 
site is the case for low n. The parameter k is linked to n to en- 
sure that half of the light is enclosed within the effective radius. 
The software GALFIT uses a 2D Sersic profile allowing for ellip- 
tical isophotes. The output from fitting the Sersic profile is given 
in terms of the total flux from the integrated profile, the effective 
semi-major axis, a e , the Sersic index, n, the axis ratio of semi- 
major and -minor axes, b/a, and the angle of a e with respect to 
the image-axes. We fitted the galaxy allowing all parameters to 
vary freely. This resulted in the following best fitting parameters: 
magAB = 24.29 ± 0.04, n = 0.95 ± 0.15, a e = 5.80 ± 0.30 px, 
b/a = 0.17 ± 0.02, and PA = -26.45 ± 1.34° in the F606W 
image. 

Due to the broader PSF in the NIR images the QSO light is 
spatially overlapping with the galaxy. We therefore fixed the struc- 
tural parameters in the Sersic fit of the galaxy to those of the well- 
constrained F606W fit. We note that simply locking all structural 
parameters of the Sersic profile in the NIR images might not give 
the most accurate description of the galaxy, as parameters such as 
size and n depend on wavelength and thus yield different results 
when analysed in different wavelength band-passes ( [Kelvin et al.| 
|2012| . However, in order to obtain the most reliable fluxes and to 
get the fit to converge we had to keep the variables fixed. In order 
to estimate how robust our obtained fluxes are with respect to the 
parameters that were held fixed, we varied the Sersic parameters 
within their errors (as given by the fit to F606W) and re-did the fit 
for each new set of profile parameters. The uncertainty on the flux 
was very minor (0.03 dex) compared to the large uncertainty caused 
by the PSF subtraction (^0.2 dex). All the obtained magnitudes are 
listed in Table |3] 

In Fig. [6] we show the WFC3 images in the three filters and 
the residuals after subtracting the modelled QSO PSF. All images 
are rotated to have North up and East left. The galaxy causing the 
DLA seen in the QSO spectrum is located to the North-East of the 
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Table 3. Results from Galfit analysis 



FILTER 


mag [ab] 


a e 1 kpc 


n 


b/a 


F606W 


24.29 ± 0.04 


1.12 ±0.06 


0.95 ±0.15 


0.17 ±0.02 


F105W 


24.51 ±0.21 


1.12* 


0.95* 


0.17* 


F160W 


23.53 ±0.13 


1.12* 


0.95* 


0.17* 



Parameters that were fixed in the corresponding fit. 



F606W 


4 



F105W 




m 



F160W 





* ♦" 



Figure 6. Data images (left) and PSF-subtracted residuals (right) from 
Galfit in all three filters. All images are aligned North up and East left. 
The galaxy causing the absorption is clearly visible in the F606W image 
even before PSF subtraction. The same position is shown by the red circle 
in the IR images. 



QSO. The red circle in the F105W and F160W images indicates 
the position of the galaxy from the F606W image. The identifi- 
cation is based on the fact that the position is consistent with the 
measured offset between the QSO trace and the emission lines in 
the X- shooter spectrum. The impact parameter of the galaxy with 
respect to the background quasar is ; /74 corresponding to a pro- 
jected distance of 6.3 kpc at z = 2.354, and the angle between the 
major axis of the galaxy and the line connecting the QSO to the 
central region of the galaxy is 71 ± 2°. 



3.7 Lya emission modeling 

We observe a typical double-peaked Lya emission line, with a 
strong component redwards of line centre and a less prominent blue 
component (Fig. [3]). The large difference between the two peaks in- 
dicates that the hya photons escape through an expanding medium, 
and the purpose of this section is to investigate what constraints we 
can put on the outflow velocity, Vout • Qualitatively we can see that 
the red wing extends to several hundreds of km s _1 . However, in 
order to get a more accurate description of the emission profile we 
have constructed a semi-realistic model of the system and run Lya 
radiative transfer (RT) through it, to find the best-fitting spectrum. 
The RT is conducted using the code MoCaLaTA (Laursen et al. 
2009a b), while the galaxy model is similar to the one described in 
( [Laursen et al.|2013) . In short, the galaxy is modelled as a sphere 
of multiphase gas, with warm, neutral clouds floating in a hot, pre- 
dominantly ionised medium. A similar procedure was used in |No-| 
terda eme et al.|p012) , although here we employ a more rigorous 
approach. 

A high number of parameters dictate the outcome of such a 
simulation. Luckily, our observations offer excellent constraints on 
many of these: The effective radius, r e — 1.12 kpc, used to model 
the size of the emitting region, the metallicity Z ~ O.31Z (we 
assume that the amount of metals that condense to dust is simi- 
lar to the local Universe (Zafar & Watson 2013)), and an average 
of 130 km s _1 from the measured emission-line widths are used 
as a proxy for the intrinsic Lya line width. From the widths of 
the low-ionisation absorption lines, we use a velocity dispersion 
of the clouds of 115 km s _1 . Finally, we infer an intrinsic equiv- 
alent width of 150 A, in accordance with the F606W magnitude. 
For details on the rest of the parameters, e.g., cloud size distribu- 
tion, temperatures and densities of the two phases, etc., see Laursen 
|eTaLl ( [20T3l >. 

We set r ga i = 10 kpc, but note that the exact size of the sys- 
tem is not important, rather the total column density (Nm) from the 
centre and out, averaged over all directions, determines the shape of 
the spectrum. This leaves us with two unknown parameters, (Nm) 
and Vout, where (Nm) is dominated by the number of clouds (N c \ ). 
We first run a rough fit to the spectrum, providing us with informa- 
tion about the initial conditions for the system: N c \ = 10 5 and 
Vout = 150 km s _1 , and consequently run a grid of 1 1 x 1 1 models 
with N c i e [10 4 - 5 , 10 5 - 5 ] and V out 6 [100, 200] km s" 1 . 

Instead of doing a regular \ 2 minimisation, we compare the 
observed and the simulated spectra using the following four ob- 
servables: The peak separation in A; the width of the red peak; 
the ratio of the integrated flux in the two peaks; and the ratio of 
the peak heights. The best-fitting outflow velocity is found to be 
Vout = 160^4q km s _1 . Here, the errors are given by the values of 
all models, for which all four estimators lie within the 68% confi- 
dence interval of those observed in our spectrum. In the best-fitting 
model, the average column density and number of clouds inter- 
cepted by a sightline towards the QSO is log(iVki) = 20.23iJ;?5 
and 2 ± 1, respectively. This is somewhat lower although not incon- 
sistent with the measured value of 20.65 db 0.05 and the fact that 
five absorption profiles were used in Sect. |3.5| 

The escape fraction of Lya photons (/ esc ~ 90%) is higher 
than what is found when comparing the total Lya-to-Ha ratio, and 
we find that a SFR of only 6.0 M yr _1 is needed to match the 
observed spectrum. However, the galaxy was modelled as a spher- 
ically symmetric system, and the galaxy appears disc-like from the 
Galfit analysis in Sec.|3.6| For such a disc-like system, the escape 
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fraction will be significantly lower when observed edge-on, which 
indeed seems to be the case here. 



3.8 Star Formation Rate 

The F606W flux of the galaxy as measured from the HST data cor- 
responds to a rest frame wavelength of 1775 A at the redshift of 
the galaxy. Based on the AB magnitude obtained from our Gal- 
fit analysis (see Table [3) we get a luminosity after correcting for 
Galactic extinction. We use the emission redshift, z em = 2.3537, 
to compute the luminosity distance. This gives a luminosity of 
L v = F v And\ (1 + z)- 1 = 9.63 x 10 28 erg s^Hz" 1 . We 
use the extinction correction factor for the F606W band, Amag = 
-0.11 mag, from the NASA/IPAC Extragalactic Databas^J(NED). 
This luminosity corresponds to a star formation rate of SFRuv = 
13.5 ± 0.5 M yr -1 using the relation from |Kennicutt| fl998). 

For the Ha emission line we take the observed line flux (see 
table [TJ, which corresponds to a luminosity of Lh« = 2.40 ± 
0.10 x 10 42 erg s -1 . Converting the luminosity into SFR again 
using |Kennicutt[jl998| ) gives SFR Hc , = 18.9 ± 0.8 M yr" 1 . 

The discrepancy between the two inferred star formation rates 
indicates some degree of dust extinction. To quantify the amount of 
extinction we assume that the two measures should yield the same 
value, when corrected for dust reddening. The correction factor to 
the SFR can be expressed as: 



SFRint — SFRobs • 10 



0A-E(B-V)k(\) 



Requiring that the two SFRs from UV and Ha be equal we arrive 
at the following expression: 



E(B -V) = 



2.5 / SFRhcA 

fc(UV) - A; (Ha) * ° g V SFRuv/ ' 

where fc(UV) and A; (Ha) denote the |Calzetti et aL^PPO) extinc- 
tion curve evaluated at the UV rest frame wavelength, 1755 A, 
and at the rest frame wavelength of Ha, respectively. From this 
relation we derive a colour excess of E(B — V) — 0.06 ±0.01 
consistent with the previously mentioned measure from the B aimer 
decrement (sect. |3.3] >. The extinction corrected star formation rate 
is SFR = 22.8 ± 1.2 M yr -1 with an assumed Salpeter initial 
mass function (IMF). If we convert this to the Chabrier IMF we get 
SFR c h = 12.7±O.7M yr -1 . This measurement agrees very well 
with the results of |Fynbo etal.| ( |2010| > and |Peroux et al.| ( |2012"] > who 
find SFRs of SFR > 10 M yr^and SFR = 17.1 ± 6.0 M© yr _1 , 
respectively. 

We note that the two proxies for star formation do not trace 
star formation on the same physical time-scales, and therefore do 
not necessarily have to yield the same measured quantity. The Ha- 
line primarily traces the ongoing star formation responsible for ion- 
ising the H Il-regions, whereas the UV continuum traces previous 
star formation as well, linked to the O and B stars already in place. 
However, the young, UV-bright stars may also contribute to pho- 
toionisation of the gas, and the quantification of the difference de- 
pends on assumptions about IMF, star formation history, and the 
distribution of dust. The detailed, exact modelling of all these fac- 
tors is beyond the scope of our simple extinction estimate, we have 
therefore chosen to neglect this effect in the above analysis. 

The star formation rate can also be determined directly from 
the absorption lines associated with the DLA if the C II* absorption 



http://ned.ipac.caltech.edu/ 



line is available (the method is described in Wolfe et al. 2003 ). Un- 
fortunately the line is blended, so we cannot put any firm constraint 
on the star formation using this method. 



3.9 Broad-Band SED Fitting 

We fit the three broad-band photometric points from the HST imag- 
ing mentioned in Sect. |3.6| to obtain estimates of stellar mass, age 
and star formation rate. But first, we apply a correction for Galactic 
extinction available from the NED online database 2 . To simplify 
the fit we subtract the emission-line fluxes from the broad-band 
fluxes and fit the continuum only, instead of fitting both continuum 
and emission lines simultaneously. The only filter which is influ- 
enced by strong emission lines is the F160W band, which contains 
flux from Hf3 and the two [O III] lines. In order to subtract the emis- 
sion lines, we assume a flat continuum spectrum, since we only de- 
tect the emission lines, and subtract the integrated line fluxes from 
the total flux in the observed band, weighted by the filter curve. We 
find that the emission lines contribute 33% of the total flux, and we 
infer a corrected F160W magnitude of ttiiqq — 23.94 ± 0.19. 

Our fitting code uses the stellar population templates from 
Bruzual & Chariot (2003 ) convolved with a large Monte Carlo li- 
brary of star formation histories (exponential plus random bursts) 
assuming a Chabrier (200 3]) IMF. Dus t is added following the two- 
component model of Chariot & Fall (2000), with the parameters 
being the total optical depth, r v , and the fraction of dust contributed 
by the ISM, yijj We restrict the range of metallicities of the mod- 
els to be consistent with our measurement from the emission lines 
log(Z/Z©) = — 0.3 ± 0.1. We then adopt a Bayesian approach 
by comparing the observed photometry to the one predicted by all 
the models in the library, and we construct the probability density 
functions of stellar mass, mean light-weighted stellar age, and star 
formation rate. Taking the mean and 16 th and 84 th percentiles of 



the PDF we obtain a stellar mass of M* 



2.1; 



x 10 y M , 



age of the galaxy of t 



98Z 



Myr, and a star formation rate 



of SFR = 8.4 +i4 M yr -1 . We also get an estimate of the dust 
extinction from the fit by comparing the intrinsic template with the 
best fitting reddened template. From this we infer Ay = 0.08+o'.o7» 
which corresponds to E(B -V) = 0.02±g;[g. 

The SFR from the SED fit agrees well with the one inferred 
in Sect. 1 3. 8 1 within la, and also the median value of the dust is 
consistent with what we find in Sect. |3.3| and |3.8| The photometry 
and best fitting template are shown in Fig. [7] 



4 DISCUSSION 

4.1 Abundances 

We have presented a measurement of the gas-phase metallicity in 
the H II regions of the emission counterpart of the DLA towards 
the quasar Q 2222— 0946. The metallicity we infer using the R23 
diagnostic expressed in Solar units is [O/H] = —0.30 ± 0.13. 
From the absorbing gas located 6.3 kpc away, |Fynbo et al. | (20 1 ) 
find a metallicity in the neutral ISM from sulfur of [S/H] = 
— 0.6 =b 0.1, thus slightly lower than the metallicity inferred from 
the central emitting region. Sulfur is not depleted onto dust, and 



For details on the prior distribution of the SFH and dust parameters see 



ISalimetal. (2005 ) 



Stellar population of a z = 2.35 DLA galaxy 9 




Viyiww?^ 



10000 12000 14000 

Wavelength / A 



16000 18000 



Figure 7. Best fitting template (in grey) to our three broad-band photomet- 
ric points (black, squares). The error-bars in y-direction indicate the \o- 
uncertainty on the fluxes. In the bottom of the figure the transmission curves 
for each filter are shown, and the blue, round points indicate the model mag- 
nitudes calculated in each filter passband. The median age and mass of the 
fit is t = 98 Myr and M* = 2.1 ■ 10 9 M , respectively. 



thereby traces the overall metallicity very well. We find a con- 
sistent metallicity from the new analysis based on our deep X- 
shooter spectroscopy. Here we found [Zn/H] = — 0.51 ± 0.06 and 
[Si/H] = —0.53 ± 0.05. Neither Zn nor Si are significantly de- 
pleted onto dust and hence provide a good measurement of the gas- 
phase metallicity ( [Meyer & Roth|1990[|Pettini et al.|1997) . From 
the elements Fe, Ni, and Mn we clearly see that dust depletion is in 
fact at play, their metallicities are [Fe/H] = -1.0, [Ni/H] = -0.9, 
and [Mn/H] = —1.2, indicating that the refractory elements are to 
some degree removed from the gas-phase, which is to be expected 
in a high-metallicity system as this particular one ( [Ledoux et al.| 
^003) . 



Bow en et al.| (|2005) find similar results regarding gas phase 



metallicities for a low-redshift galaxy, where a quasar intersects the 
galaxy 3 kpc from a star forming region. The authors find consistent 
metallicities based on the emission region and the absorbing, neu- 
tral gas, respectively. For a compilation of systems with emission 
and absorption based metallicities see Per oux et"aL| ( |2012| ). 

The presence of highly enriched material (close to Solar 
metallicity) 6 kpc above (almost perpendicular to) the galactic 
plane of this galaxy with nearly the same metallicity as the star 
forming regions within the galaxy indicates that metal-rich mate- 
rial has been expelled from the galaxy into the halo. We see inde- 
pendent evidence for outflowing gas from the Lya emission line 
with a velocity of 160 km s _1 , see Sect. 3.7 for details. At this 
velocity we estimate that it would take of the order of 40 Myr for 
this enriched material to reach a distance of 6.3kpc from the galaxy 
plane. Given this relatively short time-scale it seems reasonable that 
the two metallicities are similar, since no large amount of enrich- 
ment has had time to occur after the expulsion of the outflowing 
gas, which is expected to mix with lower metallicity gas further 
out, lowering the observed metallicity in the absorbing gas. This is 
exactly what we observe. 

A recent study by |Bregman et al.| (|2013) shows evidence for 
a similar scenario, where near-Solar metallicity gas of an edge-on 
disc galaxy has been detected 5 kpc above the disc, suggesting a 
so-called galactic fountain scenario. 



4.2 Dynamical Mass 

We can use our information about the size of the galaxy and the 
kinematics, as probed by the emission lines (see FWHM measures 
in Table [TJ, to get an estimate of the dynamical mass of the sys- 
tem. We follow the method described in Rhoads et al. (2013) to 
estimate the dynamical mass given the measured size and velocity 
dispersion: 



M dy 



<3sin 2 (i) 



where i denotes the inclination of the system with i = 90° 
being edge on. In order to estimate the velocity dispersion of the 
system we use the FWHM of the emission lines as a probe of the 
integrated gas-kinematics of the system. We then take a weighted 
mean of all the measured line-widths, and correct the FWHM for 
the instrumental resolution (45 km s _1 ). This gives a measure the 
velocity dispersion of a = 49.1 ± 7.7 km s _1 , and we get the size 
from the GALFIT analysis in physical units: a e = 1.12 zb 0.06 kpc 
(see Table [3}. 

From the GALFIT analysis we infer an axis ratio of the galaxy 
of b/a = 0.17. The system may be described as reasonably disc- 
like, given the elongated shape, and the fact that we see a value of 
Sersic n very close to 1 . Using the results of Haynes & Giovanelli 
(1984) that disc galaxies on average have axis ratios around 0.1 — 
0.2, we conclude that the galaxy in our study is very close to edge 
on, even when assuming the lowest intrinsic value of 0.1. We thus 
adopt a value of i = 90° and use the fitted half-light semi-major 
axis for our estimate of the dynamical mass of the system: Md yn ~ 
2.5 x 10 9 M©. This estimate should only be considered a rough 
approximation (valid within a factor of ~ 2) as we have assumed 
the system to be in virial equilibrium, which it may not be. 



4.3 Stellar Mass 

From our SED fit to the broad band imaging data, we obtain a stellar 
mass of M* = 2.1 x 10 9 M . We can use this measurement to test 
the recently proposed evolving mass-metallicity relation for DLA 
systems ([Ledoux et al.|2006| see also |M0ller et al.|(|2004 



|Neeleman et al. ( 2013| >). Using the relation in |M0ller etld. 



2013 



2013 



which relates emission metallicity and redshift to stellar mass, with 
our direct measurement of the emission metallicity, we find a stellar 
mass of M* = 6 1 9 x 10 9 M . Though the scatter in their relation 
is substantial (~ 0.38 dex), the agreement between the relation and 
our best fit stellar mass from the SED fit is striking. 

Moreover, the mass is in very good agreement with the median 
stellar mass of Lyman Break Galaxie s (LBGs) at redshift z = 1 — 3 
|Erb et al.|2006| . 



Hathi et al. 



( |2013| ) find M*, LBG = 2 x 10 9 



M, 



and further characterise the general LBG population in terms of 
median age (~125 Myr), median SFR (~15 M yr _1 ), and me- 
dian dust extinction ( E(B — V) w 0.15). All, but the median dust 
extinction, agree well with our inferred quantities. The dust extinc- 
tion value is, however, still consistent with our value within the 
errors (see figure 6 of Ha thi et al.|2013) . The consistency between 
the galaxy in our study and the general LBG population shows that 
there is indeed some overlap between galaxies causing DLAs and 
star-forming galaxies selected with the Lyman break technique at 
redshifts z = 2 - 3 ( |Fynbo et al.|2003l|Rauch et al.|2008) . Never- 
theless, it is important to remember that the DLA in question was 
specifically chosen to have high-metallicity, and is of unusually 
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high metallicity; it is an order of magnitude higher than the me- 
dian metallicity of DLAs at z = 2 - 3, [M/H] ~ - 1.5 or 1/30 of 
Solar ( [Noterdaeme et al.|2008[|Rafelski et al.|2012). The pro posed 
mass-metallicity relation for DLA galaxies (Ledo ux~et al.|[2 006), 
then supports the original suggestion by Fynbo, M0ller, & Warren 
( 1999 ) that most DLA galaxy counterparts are too faint to be iden- 
tified via their stellar or nebular emission. 



(2009 ) find that LBGs at z « 2 are offset by 0.4 dex compared 
to the local TFR, and they find that the relation has low scatter. 
Gnerucci et al. (201 1 ) find similar results in their sample of z ~ 3 
galaxies with an offset up to 1 dex, however, the scatter in their 
data is very large. This may indicate that the TFR has not yet been 
established at these redshifts, and that the galaxies are influenced 
heavily by random motions. 



The Tully-Fisher Relation 

We now turn to look at how this galaxy is located on the stellar- 
mass Tully-Fisher relation (M*-TFR) to test our assumption that 
the system is disc-like and relaxed. The Tully-Fisher relation, orig- 
inally stated in terms of luminosity and velocity (Tull y & Fisher| 
|1977| >, can also be presented in terms of stellar mass (which cor- 
relates with luminosity) and velocity. The M*-TFR was studied in 
detail by Kass in et aL[f2007| who gave their best fit to the data as: 



log(5 .5) = 1-89 ± 0.03 + 0.34 ± 0.05 • log 



M* 



10 10 Mr 



KV r 2 ot + <rl 



with 



where So. 5 is defined by the authors as Sk 
K = 0.5. Using the measured velocity dispersion as a proxy for 
5*0.5 (see a discussion of this in Rhoads et al. 2013 ) we find that the 
inferred stellar mass is M* = 2.6^2 x 10 9 M©. This is in very 
good agreement with our previously mentioned mass estimates, 
including our rough estimate of the dynamical mass. 

The stellar (and dynamical) mass inferred from the emission 
line widths is subject to uncertainties caused by the fact that we do 
not know the detailed structure of the velocity field. The emission 
lines are most certainly influenced by turbulence in the gas, which 
we cannot quantify. This would overestimate our line-widths and 
thereby our TF-based stellar mass and dynamical mass estimates. 
Also, we required that the system be in virial equilibrium; however, 
we observe gas at large galactocentric radii, almost perpendicular to 
the disc, with similar metallicity as the line emitting region, indicat- 
ing outflowing gas from the central parts. This has an impact on our 
calculation of the dynamical mass, and our ability to use the line- 
widths as a tracer of the ordered rotation of the system. Moreover, 
a non-negligible gas mass is expected in a young (~ 100 Myr), star 
forming galaxy. This mass would not be accounted for in the stellar 
mass estimate from the SED fit, but would contribute to the dynam- 
ical mass, thus increasing the observed velocity dispersion. We es- 
timate the gas mass from the star formation density (see Ken nicutt] 
1998) using the half-light radius and the axis-ratio of the galaxy. 



We infer a gas mass of M g 



lln * x 10 9 Mo. This estimate is 



a very rough approximation given the large scatter (^0.3 dex) in 
the relation, and is therefore not to be trusted as a true value of the 
amount of gas. It does, however, indicate that a gas mass of roughly 
half the stellar mass is present. We note that the inferred dynami- 
cal mass is consistent with our best-fit stellar mass within 1 a, and 
a significant gas mass is therefore not required in order to recon- 
cile the two mass estimates, and within the (large) uncertainties, all 
three mass estimates agree well. The compact nature of the galaxy 
also means that the kinematics of the emission lines only probe the 
innermost region, and are therefore mostly sensitive to the stellar 
mass concentrated in the central region, and not the gas in the outer 
parts. 

Studies of the TFR at higher redshifts find that the relation is 
offset to lower stellar masses for a given velocity. We see a similar 
though not statistically significant trend in our data. |Cresci et al.| 



5 SUMMARY AND CONCLUSION 

We have presented our analysis of a high-redshift galaxy selected 
from its neutral hydrogen absorption seen in the spectrum of a 
background quasar. We have presented the extracted emission lines 
from the galaxy counterpart of the absorption, and combined these 
with our detailed absorption-line study, to probe the metallicities 
seen in the two phases. We find that the two metallicities are sim- 
ilar, but the absorbing gas has a slightly lower metallicity than the 
emitting gas. We use HST imaging to constrain the stellar popula- 
tion; our data are consistent with the picture in which the galaxy 
causing the Lya absorption, is a young, small (~ 1 kpc), disc-like 
system, with a not fully ordered (proto-) disc structure. Moreover, 
we see evidence for a so-called galactic fountain, where enriched 
gas gets blown out from the star forming regions, eventually cools 
down to form the neutral hydrogen absorption that we see ~ 6 kpc 
above the galactic plane, and in the end may settle back onto the 
disc. 

This galaxy demonstrates exactly how star-forming galaxies 
at high redshift may overlap with the population of the most metal- 
rich DLAs. However, the very faint nature of damped Lyman-a 
absorbing galaxies renders the majority of these almost impossible 
to detect. The few exceptions, such as the case reported here, thus 
offer a rare glimpse into this elusive galaxy population. 
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